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Abstract 

The solitons of the SO (3) gauged Skyrme model with no pion-mass potential 
were studied in Refs. [Q, ||]. Here, the effects of the inclusion of this potential 
are studied. In contrast with the (ungauged) Skyrme model, where the effect 
of this potential on the solitons is marginal, here it turns out to be decisive, 
resulting in very different dependence of the energy as a function of the Skyrme 
coupling constant. 
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1 Introduction 

The solitons of the SO (3) gauged Skyrme model with no pion mass potental were 
studied in detail numerically in Refs. |], ||. The model 

Ho = -\^r\F l3 \ 2 + Tr ((U^D.U) 2 + k^U^D^, U^D.U}) 2 ) (1) 

was determined by the vector gauging prescription defined by the covariant derivative 

DiU = diU+ [A U U] , U = cos/ + m.<rsin/ , (2) 

where U is an element of SU(2), and the (antihermitian) Yang-Mills connection and 
curvature are in the algebra. 

Replacing dJJ in the usual Skyrme model by (@) yields the model (|l|) studied 
in P], 0| . In the present work, we augment this system 7i by the familiar pion-massQ 
potential 

7 = A(l-cos/) (3) 

such that the system we study here is TC = TCq + A(l — cos/). 

The effects of this, as we shall see from our numerical work in Section |3|, are 
considerable, in contrast with the usual (ungauged) Skyrme model || where the only 
effect is the quantitative one of rendering the asymptotic decay exponential rather 
than power like. 

We are interested in unit Baryon number solitons, hence only in spherically sym- 
metric solutions. Imposition of spherical symmetry leads to 

A = ~ {^-^ — "J ZijkVjXk , rii = Xi f = f{r) . (4) 
The static Hamiltonian reduces to the one dimensional energy density functional 
H = A(a') 2 + 2{a \~ 1)2 +l[r 2 (f) 2 + 2a 2 sm 2 f] 



r 2 2' 



+ Ana 2 sin 2 / 



(/ ')2 + a 2 Sin2 / 



+ A(l-cos/)r 2 , (5) 



2r 2 

whose equations of motion are to be solved subject to the asymptotic conditions 

lim/(r) = 7r , lim f(r) = , (6) 



lim a(r) = 1 , lim a(r) = ±1 . (7) 

Recall that in the A = case in , the boundary values (0) are instead 

lim air) = 1 , lim air) = 0, ±1 . (8) 



1 Othcr potentials can also be employed if exponential localisation is not insisted on. In particular 
for solitons with nonvanishing magnetic flux W exponential localisation is not obtained with the 
pion-mass potential 
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The asymptotic values (^j) of the function f(r) are consistent with analyticity at 
the origin and finite energy at infinity. The latter condition is ensured by the presence 
of the pion mass potential (|3]). This contrasts with the ungauged model (resulting 
from putting a(r) = 1 in for which the asymptotic condition (??) follows from 
the finite energy condition whether or not A = 0. If in the gauged model one puts 
A = 0, then condition (|6|) can be imposed as a constraint for unit Baryon charge, so 
that the solutions studied in |l], || are constrained solutions.. 

Concerning the asymptotics of the function a(r), the first member of (0) ensures 
differentiability at the origin, but the second member does not follow from the re- 
quirement of finite energy by a naive inspection of the functional (|5|). Doing the latter 
would yield instead the weaker condition 

lim a = =^> lim a = ao G R (9) 

1 — >oo r— »oo 

fixing only the r derivative of a[r) at infinity. The second members of (|7|) and (f|), 
result from a careful asymptotic analysis of the second order equations of motion 
to be given in Section [|. The results of the numerical analysis will be reported in 
Section ^ and summarised in Section [4| 

2 Asymptotic analysis: dominant balance 

The Euler-Lagrange equations for (|5[) are 

(r 2 + 8na 2 sin 2 /)/" + 2rf + lQnaa sin 2 //' + 8na 2 sin / cos f{ff 



a 4 sin 3 / cos / 



-2a 2 sin / cos /- 8k -| --Ar 2 sin/ = 0,(10) 

a" + 0(1 ~ ° 2) - ^asin 2 / - ^asin 2 f(ff - K £*£l = 0,(11) 
and the following boundary conditions are imposed: 

/—*•(), a -> a G R, /' -> 0, a' -> as r -> oo . (12) 

We want to show that there are no solutions to the differential equations (|TDD- (|TT|) 
in the asymptotic region if ao = ±1 is not satisfied. The A ^ case of interest here 
differs from the case where A = 0. First, here we will find exponential rather than 
power decay of the function /. Second, the asymptotic behaviour not ruled out by 
the asymptotic analysis presented here is ao = ±1, whereas in the case A = it is 
ao = 0, ±1. The numerical results are, of course, in accordance with the asymptotic 
analysis. 

To identify the dominant terms in the asymptotic region we use that 
8/ca 2 sin 2 //" < r 2 /", 16reaa' sin 2 //' < 2rf, 



gij^^ f cos f 

8na 2 sin / cos f(f') 2 < 2rf, 8k < 2a 2 sin / cos / < Ar 2 sin / . 
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Therefore to leading order (01]) reduces to 

r 2 f" + 2rf - Xr 2 f = 



(13) 



which can be solved in terms of Bessel functions. The asymptotic behaviour of the 
solution is 

^exp(-v / Ar) . (14) 
r 



/ ~ — exp(— y/Xr) . 

r 

Note that in the case A = 0, equating the terms of leading order in ( fTO|) yields 



r 2 f" + 2rf-2a*f = 0. 



This equation has 



/ = / or -(l+V 1+8a o)/2 



(15) 
(16) 



as a family of solutions which satisfies the boundary condition for /. 
In the asymptotic region we also have 



1 



sin /(/') <C -a sin /, k 



'3 • 4 £ 

a sin j 



<C -a sin 2 / . 
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In terms of A = a — ao we therefore obtain from flTT| ) to leading order 



r 2 A" + (1 - 3a 2 )A = a (a 2 - 1) + — (a + A)f . 



;i7) 



Because of the exponential decay of / in the case A ^ 0, the equation reduces to 

r 2 A" + (1 - 3ag)A = a (ajj - 1) . (18) 



For 3<2q = 1 the general solution of (JlS 



is 



.4 



3^ 



logr + Ci + C2T 



(19) 



which does not go to zero and is unacceptable. The solutions of (|T8| ) for \a Q \ < | 
a = ±i and \a \ > | are, respectively, 



,4 



a — ao 
1 -3a 2 , 



+ Vr 



c 3 cos(- ^3 - 12a 2 , logr) + c 4 sin(-y / 3 - 12a 2 , logr) 



A = =f- + c 5 y/r + c 6 \/r\ogr 



A 



ftp ~ clq 
1 — 3a 2 , 



+ \Jr c 7 r 



,^12^-3/2 



+ c 8 r 



yi2ag-3/2 



(20) 
(21) 

(22) 



which for ao 7^ ±1 do not go to zero either. So there is no solution of the equation 
flnp for which ao 7^ ±1. Equation fllBD has solutions with acceptable asymptotic 
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behaviour only if a = ±1, given by (P2|) . In this case A ph Cgr" 1 (with cj = 0)is a 
solution with acceptable behaviour in the asymptotic region. 

The situation here should be contrasted with the case where A = 0. There [I], for 
a = 0, 

/«^, A « Aor^v 7 ^)/ 2 (23) 
r 

is an acceptable solution of ([T3|) and ([H]) in the asymptotic region if |/o| > 2. In all 
cases the existence of acceptable solutions for small and large r of course does not 
guarantee the existence of a solution which has both, acceptable behaviour for small 
and for large r. 

3 Numerical results 

Before presenting our numerical results for the A ^ case, let us recall those for the 
system (|5|) with A = 0, carried out in [[I], |2|. With respect to the Skyrme coupling 
constant k, only two energy branches of stable solutions were found. These are 
presented in Figure 1 and labelled A and B respectively. The branch A corresponding 
to a(oo) = 1 starts from k = and persits up to a cusp critical point . The classical 
energy is an increasing function of k (with E(k = 0) =0). Branch B corresponding to 
a(oo) = starts from another critical point K^ r and persists up to k — > oo. Branches 
A and B are bridged by a branch A' which like A has a(oo) = 1 and is known to be 
unstable 0; the numerical values nf ~ 0.8091 , k c { ~ 0.6914 were found in 0. No 
solutions with a(oo) = — 1 were found, for which we do not have an explanation. 

As we have no explanation for our inability to construct solutions with a(oo) = — 1 
for the A = system, we have used a different numerical procedure as well as the 
alternative ways to implement the boundary conditions to attempt to construct them 
and/or to recover the results of Jj], 0. Namely, this involves the use of the boundary 
values (U) for the chiral function f(r), while for the gauge function a(r) we have 
used the first member of (0) together with (|9|). In this way, we indeed recovered the 
branches A, A' and B found in [p], |2| but nothing else. This provides a good check 
on the correctness of the numerical procedure to be applied to the A ^ system in 
the present work. 

After all these tests involving the boundary condition (^) we have applied our 
numerical routines to the A ^ system and we now discuss the results. Various types 
of solutions are presented below for the value of A = 1.0 but we have checked that 
for generic values of this parameter, the corresponding solutions display the same 
qualitative features . 

First, the numerical analysis confirms the occurence of a branch of solutions obey- 
ing a(0) = a(oo) = 1. It exists for all values of k and the energy increases mono- 
tonically with this parameter, as illustrated on Fig. 2 (branch ^4). This solution can 
further be characterized by x m , the value of x where a(x) attains its local minimum. 
The following results hold : 

lim K ^o£m = , Um^ooXm = x m w 1.33. (24) 
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For sufficiently high values of k, our numerical analysis also reveals the existence of 
two branches of solutions obeying a(0) = — a(oo) = 1. For instance, for k > k c 
(k c ~ 2.297 in the case A = 1), we constructed two different branches of solutions, 
labelled C an C on Fig. 2. As indicated the classical energies of these solutions are 
slightly higher than the corresponding one on branch A. In the limit k —>■ k c , the two 
branches C and C coincide. 

We also find (using obvious notations = energy on branch A,...) 

Mm^Ec - E A ) = . (25) 

The two new solutions can further be characterized by the value x$ where the function 
a(x) takes its (unique, as far as we can see) node. The evolution of function 
of k for the two branches is displayed on Fig. 3; this figure clearly indicates that, 
for k fixed, the two solutions are distinguished by the position of their node. The 
numerical results suggest that 

lim^ooZo = lwi K -+ 00 x m = x m for branch C , (26) 

lim^ooXm = oo for branch C , (27) 

Finally, on Fig. 4, the profiles of a(x) for different values of k and on the different 
branches are compared. It is seen that, on the interval [0, x m ], the function a(x) 
corresponding to the solutions C and C deviates only a little from the corresponding 
one on branch A . The differences between the three solutions are rather perceptible 
on the interval [x m , oo]. 

Quantitatively the value of A affects the value of k c slightly. We find numerically 
k c = 2.423 , 2.359 , 2.297 respectively, for A = 0.1 , 0.5 , 1.00. For all cases, we 
checked that the profiles of a(x) and of f(x) obey the asymptotic behaviours obtained 
in the previous section, providing a nice check of our numerics. The various solutions 
were obtained by using the subroutine COLSYS || (based on the damped Newton 
method of quasi-linearization, a brief description of it is presented in the Appendix 
of 0). The solutions were obtained with a high degree of accuracy : typically with 
errors less than 10 -8 . 

The solutions displayed by Fig. 2, constitute the main result of the present work. 
Comparison of Figs. 2 and 1 shows that the inclusion of the pion-mass potential in 
the gauged sigma model has a considerable qualitative effect compared to the case 
when this potential is absent. 

Concerning the stability of the pattern of solutions displayed in Fig. 2, we have not 
carried out a detailed quantitative analysis. However we beleive that the solutions on 
the branches A, C and C possess respectively zero, one and two unstable modes. The 
stability of branch A is guaranteed by the toplogical lower bound associated with the 
nontrivial asymptotics of the Skyrme field. The plot of the energy corresponding to 
the branches C, C terminates into a cusp, typical of catastrophe theory (see e.g. ||). 
In such situations it is believed (and it was demonstrated numerically in a particular 
case |§) that the number of negative modes of the upper branch of the cusp exceed 
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by one unit the number of negative modes of the lower branch. With our expectation 
of the number of unstable modes, the calculation of the Morse index 

£ = £("! W (28) 
q 

(the sum runs over the classical solutions, q counts the number of negative modes of 
the solution and N q represents the number of them) leads to £ = 1 irrespectively of 
the parameter k. 

4 Summary 

We have investigated the effect of adding the pion mass-potential, namely the last 
term in @, to the Skyrme model. This potential enforces the conventional asymptotic 
values (H) of the chiral function /(r), but in the gauge- decoupled case, with a(r) = 1 
dD) (i) i s independently guaranteed by the finite energy condition. Not surprisingly in 
that case ||, with a(r) = 1, switching A on or off results in no appreciable qualitative 
change in the Skyrme soliton. The power dacay for the A = case is replaced by an 
exponential decay in the A 7^ case. 

The situation is quite different in the gauged case, where the potential is necessary 
to have (|6p. On the other hand imposition of tffi) also guarantees unit Baryon charge, 
hence it is tantamount to insisting on unit Baryon charge even if A = in (||). This 
is because the Baryon charge density g top , is a topological charge density and hence 
essentially a total divergence. Thus if one sought conditional solutions restricted to 
unit Baryon number by adding the density ^Qtop to the static Hamiltonian, with £ 
the Lagrange multiplier, then the Euler-Lagrange equations would remain unchanged, 
justifying the integration of these subject to @ without it being necessary to have 
A 7^ in (Q). The A = case was the problem analysed in 0], and the solutions 
constructed there must be considered to be conditional solutions. Here we have 
studied the case woth A ^ 0, and found that switching A on results in considerable 
qualitative effects on the (gauged) Skyrme soliton. 

This is not a surprising result, but nor is it predictable. In the present case A 7^ 
with asymptotics (|^), it was possible to construct solutions numerically with both 
asymptotic values a(oo) = ±1. By contrast in the A = case with asymptotics (||), it 
turned out [|I|, |2| that out of the three possible asymptotic values (||) only for the two, 
a(oo) = and 1, could solutions be constructed numerically. Moreover, the respective 
energy profiles of these solutions are quite different as seen from Figs. 1 and 2. 

Concerning the phyical relevance of these solitons, first we note thay they describe 
the low energy properties of the Nucleons and then note that the SO(3) gauging 
of this system is a perfectly natural step in the direction of the studying the effects 
of the 50(3) x 50(2) Standard Model on the Nucleons. 
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Figure Captions 



Figure 1 Energy versus k for a(oo) = 1 solutions (branches A, A') and for a(oo) = 
solutions (branch B), for (for A = 0.0). 

Figure 2 Energy versus k for the a(oo) = 1 solution (branch A) and for the a(oo) = — 1 
solutions (branches C,C) (for A = 1.0). 

Figure 3 The evolution of the position of the node of the function a(x) as a function of 
k for the branches C and C . 

Figure 4 The profiles of a(x) for the three solutions available for k = 2.5 are presented : 
C in solid line, C in short-dashed line, A in dot-dashed line. The profiles of a(x) 
at the critical value (k ~ 2.3) and for k = 2.5 on branch C are supplemented 
respectively on the long-dashed and dotted lines. 
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